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ABSTRACT 

The massive black hole at the Galactic center Sgr A* is surrounded by a cluster 
of stars orbiting around it. Light from these stars is bent by the gravitational field 
of the black hole, giving rise to several phenomena: astrometric displacement of the 
primary image, the creation of a secondary image that may shift the centroid of Sgr 
A*, magnification effects on both images. The near-to-come second generation VLTI 
instrument GRAVITY will perform observations in the Near Infrared of the Galactic 
Center at unprecedented resolution, opening the possibility of observing such effects. 
Here we investigate the observability limits for GRAVITY of gravitational lensing effects 
on the S-stars in the parameter space [-Dls ; 7j K], where -Dls is the distance between 
the lens and the source, 7 is the alignment angle of the source, and K is the source 
apparent magnitude in the -fT-band. The easiest effect to be observed in the next years 
is the astrometric displacement of primary images. In particular the shift of the star 
S17 from its Keplerian orbit will be detected as soon as GRAVITY becomes operative. 
For exceptional configurations it will be possible to detect effects related to the spin of 
the black hole or Post-Newtonian orders in the deflection. 

Subject headings: Gravitational lensing — Black hole physics — Galaxy: center - 
astrometry — infreared: stars — instrumentation: interferometers 

1. Introduction 

Understanding the nature and the physics of a black hole has always fascinated generations of 
scientists as well as ordinary people. Obviously, the first basic question is whether there exist black 
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holes actually. Only recently the astronomers have been able to gain quite conclusive evidence to 
support a positive answer. For instance, it is widely accepted that the central engine of an AGN is 
an accreting supermassive black hole (M. > 10 6 Mq) at the center of the galaxy, responsible of the 
observed high-energy emission; stellar black holes (3M < M» < 10 Mq) are observable in close 
binary systems, when matter is transferred from a companion star to the black hole; many ultra- 
luminous X-ray sources are finally suspected to be intermediate-mass black holes (M. > 10 Mq). 
However, the royal evidence of the black- hole existence resides in the innermost part of our Galaxy, 
and comes from decennial near-infrared (NIR) observations of the orbital motion of more than two 
dozen stars (the so-called S- stars) around the well-known and mo tionless compact radio source, Sgr 



A*, which is indeed a massive black hole of about 4.3 x 10 6 Mq (lEckart et al 



2003 



2009a 



Ghez et al 



2005a 



Schodel et al 



Eisenhauer et al. 



2005 



Ghez et al 



2008 



Trippe et al 



2002 



2008 



Schodel et al 



Gillessen et al 



201(3 ) . In particular, more than 15-years continuous monitoring of this stellar 



cluster allowed t o observe the whole or bit of the star S2, which completed a full revolution around 
Sgr A* i n 2008 ( Gillesse n et alj|2009bl ). For an exhaustive review on the argument, we refer the 



reader to Genzel et al 



(I2010h . 



Therefore, Sgr A* represents by far the best (and closest) place for searching for post-Newtonian 
effects in strong-field regime, like for instance the relativistic motions of flaring jets and infalling 
matter close to the event horizon, the deviation of stellar orbits from Newton dynamics (rela- 
tivistic periastron shift and rotation), the precession of the angular momentum vector around the 
spin of the black hole (Lens Thirring Precession), redshift variations in the spectrum of orbiting 
stars, etc. General relativity clearly represents a golden key to unravel the physical features of 
the black h ole (spin, inclination, etc.) and distinguish among the various (competitive) theoreti 



2005 



cal models (lJaroszvnskilll998bl:lFragile fc Mathewsl)200o1: blubilar &: Eckart 



Weinberg et alJl20 05: Zuc ker et al]l200fl:IWillll2008l : lBroderick et al.ll200' 



Merritt et al 



2009; 



Merritt 



20ld : lAngelil fc Sahall2011 



9: 



20011: iBroderick fcLoeb 
Hamaus et all 2009; 



Another very interesting effect expected from general relativity is gravitational lensing. In 
fact, Sgr A*, acting as a powerful gravitational lens, is able to deflect the light rays emitted by 
the S-stars from t heir trajectories, affecting; th e ir measured positi o ns and eventually generating 



Chaname et al 



secondary images (IWardle Yusuf-Zadch 



2001 



Bozza fc Mancini 



200 



5, 



1992 



200. 



3: 



Jaroszvnski 



1998a 



Bozza &: Mancini 



2004; 



Alexander &; Sternber 



Nusser &: Broadhurs 



1999; 



2004; 



De Paolis et al 

2009T). Howeve r, the expected astrometric displacements of the primary 



image (~ 20 //as, IGillessen et al.l 2009a) an d the secondary-image luminosities (K = 20.8 in the 
present known best case, iBozza &: Mancinill2009i ) are difficult to detect and the resolution of the 
most powerful modern instruments is currently insufficient to perform such high precision astrom- 
etry and photometry. However, the perspectives for observing S-star displacements and secondary 
images generated by the massive black hole at the Galactic center are not so far from what the 
modern technology is going to put into play. Innovative NIR interferometry in struments are now 
under devel opment at the Very La r ge Telescope and Ke ck: they are P RIMA (ID elplanckej 120081 ) . 



GRAVITY (jEisenhauer et al.ll2008l : IGillessen et alJl201(j ). and ASTRA (jPott et al 



2008 



These 
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instruments have been conceived to achieve an astrometric accuracy of 10 — 100 //as in combination 
with milli-arcsec angular-resolution imaging. 

Here we focus on the new second-generation Very Large Telescope Interferometer (VLTI) 
instrument, GRAVITY, that w ill be scientifically oper ative at the ESO Paranal observational site 



starting supposedly from 2014 (jEisenhauer et al.ll2008l ). Its precision will allow to study the close 
environment of Sgr A* in detail, testing general relativity in this so striking laboratory. The aim 
of the present paper is analyze the ability of GRAVITY to detect gravitational-lensing effects 
generated by the central black hole, estimating the position and the luminosity of the primary and 
secondary images of the S-stars as functions of time. 

Our paper is structured as follows. In § 2 we review the current NIR observations that have 
lead to the present knowledge of the Galactic center region. We also introduce GRAVITY and its 
expected performance. In § 3 we examine the gravitational lensing effects by dividing the space 
around Sgr A* in three regions according to the GRAVITY ability to separate the images of a 
stellar source from the black hole. In § 4, 5, and 6 we illustrate the respective results in the three 
regions. In § 7 we draw the conclusions. 



NIR observations toward the Galactic center 



From a physical point of view, the central region of the Milky Way is on e of the most interestin g 
zones of the Universe. Being located at Rq = 8.3 kpc from the Solar system (jGillessen et al.ll2009al ). 
it is the nearest and most accessible place where to study the interplay between a central massive 
black hole with the gas and the stars in its environment. This region is also characterized by a 
very high interstellar extinction (roughly 30 mag in the optical band) along our line of sight, since 
it crosses the disk and several spiral arms of the Galaxy. Nevertheless, the Galactic center is quite 
transparent in the NIR band s, where the extincti on amounts to only ~ 3 mag. In particular, a 
detailed study performed by ISchodel et al.l (|2007l ) showed that the interstellar extinction in the 
central ~ 0.5 pc is high ly variable and has a general minimum centered on Sgr A* of ^4k = 2.8 
(jEisenhauer et al.l 120051 ). 



2.1. Current NIR observations 



After the first observations by speckle imaging , which have placed the first constraints on the ex- 
istence of a mass ive black hole (jEckart et al.lll993l : iGenzel et al.lll99a : lEckart fe Genzellll99a . Il997l : 
Ghez et al.l Il998l ). the NIR observations (both imaging and spectroscopy) of the galactic center 
are currently performed thanks to the adaptive optics and laser guide techniq ue, two very innova 



five technologies introduced in the large-class telescopes (8 — 10 m) since 2002 fSchodel et al 



2002 



Ghez et alJ[20oil : ISchodel et aDbood : lEisenhauer et alJbood : bhez et al.ll2005bl . l2008l : lGillessen et al. 
2009al ). Thanks to these instrumentations the imaging observations can go about 2-3 mag deeper 
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than previous speckle imaging observations (jSchodel et al.l 120071 ). 

Considering the high density of the star cluster in the neighborhood of Sgr A*, a 8 m-class 
VLT unit telescope is able to detect a faint star down to an apparent magnitude of 18.5 in the 
K band (around A = 2.2 fim). Instead, the limit of the spectroscopical measurements allows to 
identify early- and late-type stars as faint as K ~ 16.5 and K ~ 17, respectively. For the brightest 
S-stars, it is also possible to achieve astrometric and rad ial- velocity measurements with a preci sion 
in the range 200 — 300 ^as and 15 km s -1 , respectively (jFritz et al.ll2010l ; iGillessen et al.ll201ll ). 



2.2. Sgr A* 



The center of the Milky Way hosts a very exotic source, which is variable across roughly the 
entire electrom agnetic spectrum. It h as been detected for the first time in the radio wavelength 
many years ago (|Balick Brownlll974h . and nowadays we know that this source is very compact and 
moti onless, and its location coincides with the dynamical center of the Galaxy (jReid Brunthaler 
2004). Over the last decade, Sgr A* has been fully analyzed across the electromagnetic spectrum, 
at radio, sub- millimeter, NIR, and X-ray wavelengths, revealing itself as a highly variable source 
especially in the X-ray and NIR bands. Very extensive observations in the K band taken at 
NACO/VLT in 2004 — 2009 showed that Sgr A* has a continuously variable low-level state. The 
low-level emission above 5 mJy identifies the flaring emission state and the low-level emission below 
5 mJy is considered as the quiescent state, with at least 0.5 mJy of the flux due to a faint stellar 
component. The quiescent state is variable with four NIR peaks and one X-ray peak a day. Flares 



that have more than 5 mJy of variability are quite rare, roughly two per year (|Dodds-Eden et al 

201 ih . 



2.3. S-Stars 

The central 10 arcsec region of the Milky Way hosts a stellar cluster, formed by old, late-type 
red giants, supergiants, asymptotic giant branch stars, but also by many hot young, early-type 
stars, like post-main sequence blue supergiants and Wolf-Rayet stars. In particular, deep NIR 
observations of the Galactic center (both proper-motion and Doppler measurements) revealed ~ 100 
stars within r = 1" (~ 0.04 pc) of Sgr A*, with a remarkable concentration of B stars, the so-called 
S-cluster 0. 



The existence of so many young stars in this place raises an intriguing paradox (jGhez et al 



2003). Due to the intense gravitational field generated by the central black hole, it is improb- 



able that they formed in situ, but on the other hand their age (2 — 8 Myr) is not in agree- 



Unfortunately, the MPE and the UCLA groups used different nomenclature for the stars beloging to this cluster. 
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ment with migration time scale. Other alternative scenarios have been postulated, including a 
mech anism of exchange cap t ure of young binaries and the demand of an intermediate mass black 



hole (Gould & Quillen 



2003 



Milosavlievic fc Loebll20o3: [rlarfst et al.ll2008l: lLockmann et al.ll2009l : 



Madiean et al J 12009 : iPerets et al.ll2009l : iMerritt et al.ll2009l : lGrivil2010l : lAlexanderll201ll ) 



By using the above cited high-resolution near-infrared techniques, 16 years of extensive moni- 
toring of stellar orbits around Sgr A*, have allowed to refine the orbital p arameters of many S-s tars. 
At the present time, the orbits of 27 S-stars have been well determined (Gillessen et al,ll2009ari. In 
a gra vitational-lensing context, these stars have already been analyzed bv lBozza Mancinil (|2005l . 
2009), who calculated all the properties of their secondary images, including time and magnitude 
of their luminosity peaks and their angular distances from the central black hole. Currently, the 
most interesting case came out to be S6, since its secondary image will reach K = 20.8 at its peak 
in 2062, with an angular separation of 0.3 mas from the central black hole. 



2.4. GRAVITY 



A big step in angular resolution is expected to come from NIR interferometry. Actually, 
GRAVITY is a second-generation VLTI instrument, specifically designed to observe highly rela- 
tivistic motions of matter close to the event horizon of Sgr A*. With a baseline of ~ 100 m, it 
should ensure the access to its innermost stable circular orbit (the radius of this orbit is 30 /ias 
for a 4.3 x 1O 6 M Schwarzschild black hole), allowing astrometric detection of hotspots orbiting 
around the black hole. Thus, the high sensitivity of GRAVITY should allow to reconstruct images 
with details as faint as K = 19, yielding a resolution of ~ 3 mas for objects that can be as faint as 
K = 18. With this remarkable resolution, GRAVITY should be able to catch the light of most of 
any very fast orbiting stars within the central 100 mas. Its design wi ll finish during 2012 and it w ill 
be presumably installed in Paranal in 2014. (jEisenhauer et al.ll2008l : iGillessen et al.ll2010i . l201ll ). 



According to recent numerical simulations by I Vincent et al.l (|201ll ). in the pure imaging mode, 
GRAVITY should provide an astrometric precision of the order of, or better than, the Sgr A* 
Schwarzschild radius, for a source limiting brightness of K < 15 and with an integration time of 
100 s. However, it must be taken into account that the astrometric precision is a lso very dependent 
on the number of sources in the field dGillessen et alJboiol : Ivincent et al.lboilh . 



3. Gravitational lensing phenomena around Sgr A* 

The focus of this paper is on the effects of gravitational lensing by Sgr A* on nearby stars and 
how these affect their apparent properties such as position and luminosity. We will model Sgr A* 
as a purely Schwarzschild black hole, arguing on the validity of this hypothesis for the proposed 
observations all along the illustration of the results. 
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3.1. Basics of lensing 



We assume a dista nce to Sgr A* of Dq,l = 8.3 kpc and a mass of M = 4.3 x 10 6 M , following 
Gillessen et al.l (|2009al ). We will refer to the line connecting the observer to the lens as the optical 
axis. The source position is then fixed in polar coordinates around the lens by a radial distance -Dls> 
a polar angle 7 taken from the optical axis, and an azimuthal angle eft. The latter is irrelevant in 
the lensing discussion, given the assumed spherical symmetry of the lens. 7 ranges from (perfect 
alignment of the source behind the lens) to ir (perfect anti-alignment, with the source in front of 
the lens). In Figure [U we illustrate the notations in the geometric set-up for gravitational lensing. 



The Schwarzschild radius of Sgr A* is 



Rs 



2GM 



1.27 x 10 10 m 



4.12 x l(T 7 pc. 



(1) 



The angular Schwarzschild radius is thus 

S = Rs/D OJj = 10.2 fias. 



(2) 



The motion of photons follows null geodesies around a Schwarzschild black hole. We can thus 
write an exact lens equation relating the source position (Duiij} to the image position 9 a s s een by 



the observe r (cfr. Figured]) (jFrittelli &: Newman 



2008 



1999 



Frittelli et al 



2000 



Perlick 



2004. 



3H; 



Bozza 



20101 ) . To this purpose we define 



A(r) 



Rs 



(3) 



i.e. the gu component of the Schwarzschild metric or equivalently the inverse of the g rr component. 
Considering that the impact parameter u is related to the image position by 



: s'mO ~ DqlO, 



and that the closest approach distance of the photon r m is related to this by 



(4) 



u 



A(r m y 



(5) 



the lens equation assumes the form 

2tt qp 7 : 



Sol D h § 
+ 



r^/r 2 — u 2 A(r) 



zdr, 



(6) 



where the upper sign holds for the primary image 9 + and the lower sign holds for the secondary 
image 9- appearing on the opposite side. 
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The i ntegral can b e easily evaluated in terms of elliptic functions as widely reported in the 
literature (|Darwinlll959l ). 

The magnification of a gravitational lensing image can be calculated by taking the ratio of an 
angular displacement in the observer sky and the corresponding displacement in the source surface 
orthogonal to the emission direction of the photons. Following iBozzal (|201Cj), and expressing the 
result in a synthetic form, we have 



D 2 
u OS 



D ts D OL sin 7 



1 



du du 
+ u- 



dD 



LS 



(7) 



where Dos = y D OL + D^ s + 2L>ol-Dls cos 7 is the distance from the observer to the source. The 
derivatives du/dj and du/dD^s can be evaluated numerically using the lens equation. Note that for 
small impa ct parameters co mpared to the source distance (u <C -Dls) we recover the magnification 
formula by lOhanianl (|l987l ) . In the context of black hole lensing, this is sometimes referred to as 
thin lens approximation. 

Gravitational lensing effects are maximal for 7 ~ (standard lensing) , with a second maximum 
for the secondary image occurring for 7 = ir (retro- lensing) . This is reflected by the sin 7 in the 
denominator of the magnification formula. 

In the weak field limit, the integral of Eq. @ simplifies and the lens equation becomes 



±7 



+ arcsin 



Rs 



R S (2D LS + r m ) 



Dls{D L s + 



(8) 



where we have just saved the lowest order term in r m /Z^oL; whereas we allow -Dls ^o become 
comparable to r m . 

The relation between the minimum distance r m and the impact parameter u can be written 
to first order in R$ as 

Rq 

r m = u-^+o(R s ). (9) 



In addition to the weak field approximation, the thin lens approximation requires Dls 3> r m, u. 
This automatically implies small angles, so that we can set 7 = Dos/3/ D^s, where f3 is the observed 
angular position of the source if there were no lens (see Figure [1]). The small angles weak field lens 
equation is then obtained 



where 



is the so-called Einstein angle. 



/3 = #--f 



Oe 



2R s Dls 
DqlDqs 



(10) 
(11) 
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Summing up, the full general relativistic treatment without approximation is given by Eq. ©. 
The weak field approximation is given by Eq. ([8]) and the weak field + thin lens approximation 
is given by Eq. (|10p . Note that if we want to test the sensitivity of gravitational lensing to Post- 
Newtonian orders of the black hole metric, we must compare the results of Eq. (JSj) to Eq. ([6]). In 
fact, in many cases the differences between Eq. © and Eq. (|10p are not imputable to strong field 
effects but just to the thin lens/small angles approximation. 

Nevertheless, it is useful to recall the simple analytical expressions of images and magnification 
obtained by Eq. (fTUj) . since they provide a useful guide to the expected gravitational lensing effects. 
We have 



O± = -[P±y/lP + 40*) (12) 

±11 (13) 



The primary image is always outside the Einstein ring with radius 9e, while the secondary 
image is always inside the ring. Both images approach the Einstein ring as 7 approaches 0, merging 
into a single degenerate ring image with radius 0e for 7 = 0. 



3.2. Strong field and spin corrections 

Sgr A* offers a unique chance to test GR in a strong field regime. It is thus of primary 
importance to understand to what extent a weak field Schwarzschild treatment based on Eq. (|10p 
is sufficient to explain all the phenomenology within reach of GRAVITY and at which point the 
full GR treatment of Eq. (jSJ) is really necessary. This will also indicate where to look for new tests 
of GR based on lensing observations. 



3.2.1. Black holes in alternative gravity theories 

Most alternative theories of gravity only become effective in the neighborhood of the event 
horizon, by predicting space-time metrics deviating from Schwarzschild in such an extreme regime. 
These deviations are presumed to assume sizeable proportions for Sgr A*, which is the closest 
massive black hole to us. At the same time, gravitational lensing provides one of the simplest ways 
to test these alternative theories once the metric is known. For this reason, gravitational l ensing 



by Sg r A* is invoked a s the natural astrophysical test in these theoretical studies (see e.g. iBozza 



(|2010l ): iPerlickl (|2004bl ) and references therein). This class of alternative metrics includes black 
holes with scalar fields, string theory black holes, braneworld black holes, wormholes, Born-Infeld 
gravity black hole. 
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All these metrics predict appreciable effects for the higher order images, i.e. those images 
generated by photons looping around the black hole once or more times. However, these higher 
order images require resolutions of the order of ^as or better and exceptional alignments to be 
detected, since they are extremely demagnified. In the context of lensing of S-stars, higher order 
images were examined bv lBozza Mancinil (|2004l . 120051 ) for the Schwarzschild metric, where it t is 
shown that they are always fainter than K = 30 for all known S-stars. A firs t direct attempt to fin d 



2011 



effects of alternative gravity theories in lensing of S-stars has been lead by iBin-Nunl (|2010l . 
who considered a tidal Reissner-Nordstrom metric coming from Randall- Sundrum II theories. The 
perturbations are larger for images forming very close to the black hole, but become very small for 
secondary images forming far enough from the black hole (AK ~ 0.01 for S6). 

Summing up, the common feature of all these metrics is that they reproduce the Newtonian 
limit far from the black hole, while they differ from Schwarzschild in their Post-Newtonian expan- 
sions. The position and the magnification of the images can be calculated analytically beyond the 
weak field limit of Eqs. (|12p -( [T3"l) for a generi c Post-Newtonian metric, by an expansion in powers 



of the Schwarzschild radius ( Ebina et al 



2000 



Lewis h Wang 



2001 



Sereno 2004; Keeton & Petters 



20051 ). The relative correction is of order 



56± _8fi _e s 



(14) 



It therefore takes effect only for light rays passing very close to the event horizon. This happens 
for sources very close to the black hole (e.g. on the accretion disk) or for secondary images generated 
by sources far from the optical axis. In the latter case, the secondary images also become very faint 
and difficult to distinguish from the lens. For this reason, S-stars do not seem to be quite ideal 
candidates to study strong field effects. However, we will see that there are regimes in which higher 
order deflection terms may become marginally accessible. Therefore, for each observable, we will 
indicate the region in which a source should be in order to yield sizeable Post-Newtonian effects. 



3.2.2. Alternatives to the black hole hypothesis 



Alternatives to t he standard black hole paradigm com e from boson ([Torres. Capozziello fc Lambiasd 

2000), fermion stars ( Viollier. Trautmann &: Tuppeifll993 ) and dark matter concentrations ( Tsiklauri Sz Viollier 
19981 ). While t hese objects are bui l t so a s to reproduce the observed spectrum of the black hole 
accretion disk (IGuzman &: Becerrill l2009i ) . their extended nature leaves a fundamental signature 
in gravitational lensing, since the deflection angle no longer grows when the impact parameter 
is decreased below the physical radius of these stars. However, the restrictions imposed by their 
construction scheme makes this difference arise only for light rays approaching closer than a few 
Schwarzschild radii. So, once more a source very close to the central object is needed. 
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3.2.3. Modified gravity 



Among the many attempts to solve the dark matter /dark energy problems by modifying GR, 
the MOND th e ory is the most famous (jMilgroml Il983l ) , thanks to its recent relativistic version 
by iBekensteinl ([2004), called TeVeS. In this case, gravity is modified at small fields, when the 
gravitational acceleration drops below a new fundamental c onstant an ~ 10~ 10 m/s 2 . Gravitational 
lensing by point-like objects in this context was studied by lChiu. Ko. Tianl (|200a ). It turns out 
that deviations from Schwarzschild lensing arise for impact parameters comparable or larger than 
the scale 

r = VGM/a ~ 77 pc. (15) 



These effects would thus be important for far background bulge stars lensed by Sgr A*, in 
which case the competing effect of the stellar mass of the inner clusters that would sum up with 
the black hole lensing would probably make a clear detection of TeVeS lensing very difficult. 

In the case of S-stars, whose distance from Sgr A* is below the pc, TeVes corrections are largely 
negligible (order nas at most in the astrometric shifts). 



3.2.4- Black hole rotation 



Rotating black holes are described by the Kerr metric, characterized by the spin parameter a, 
corresponding to the black hole specific angular momentum, ranging from (Schwarzschild) to M 
(extremal rotation). In Kerr black holes, the point-like caustic of the single lens becomes a small 
astroidal caustic, shifted from the optical axis by aR$/2M in the west direction (if we identify north 
as the direction of the black hole spin projected on the sky) (jSereno & De Lucall2008l ). This means, 
that in a first app roximation the Kerr b l ack hole works as a Schwarzschild black hole slightly shifted 



from its position (jAsada 8z Kasai 



2000 



Asada. Kasai Yamamotd 120031 ) . If the center of Sgr A 



is pinpointed to very high accuracy (below the angular extension of the Schwarzschild radius), then 
all gravitational lensing effects will work as if the black hole were displaced by aRg/2M. Since 
many quantities diverge as the inverse of the distance of the source to the optical axis (-Dls7 i n 
our notations), even such a small displacement of the caustic may induce dramatic effects. For this 
reason, we will also indicate the limits of sensitivity to the spin for each observable we are going to 
examine. 



3.2.5. Thin lens approximation break- out 

The last possible break-out of the weak deflection lens equation (fTQj) arises when the source 
is very badly aligned (7 3> 1). Then the primary image cannot be calculated using the so-called 
thin-lens approximation, which in this context consists in sending the integration limits to infinity 
in Eq. ([6]). In this case, in fact, Dls becomes comparable to the impact parameter u. 
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3.2.6. Image distortion 



Finally, we also note that it would be extremely improbable to see any kind of image distortion 
of the stars around Sgr A*. In fact this is possible only for alignment angles 7 of the order of the 
source angular radius, which in the case of S-stars is of the order of ns of arc. For comparison, in 
Figure [2] we report the positions of S27 and S6 at the best alignment epoch of their orbit. These 
two stars, which represent the best known cases, at most get at 7 ~ 3°. Therefore, we will always 
consider the images of the S-stars as effectively point-like for GRAVITY and thus described by a 
PSF. For the same reason, the finite size of the Kerr caustic in a rotating black hole does not lead 
to any observable effects. 



Gravitational lensing generates two images of a given source star. These images might be 
confused or not within the PSF of Sgr A* itself, depending on the angular position of the images 
compared to GRAVITY resolution limit of 3 mas. We can therefore distinguish three different 
regimes for gravitational lensing, occurring in three different regions of the plane (7, Dls)- 

• Region I: No images are resolved from Sgr A*. A single blend appears containing both images 
of the source star along with the radiation from the central black hole. 

• Region II: The primary image is resolved separately while the secondary image is still confused 
with Sgr A*. 

• Region III: Both images are distinguished from Sgr A* and can be studied separately. 

Note that \Q-\ < \9+\ so that it is not possible to have the secondary image resolved with the 
primary unresolved. 

The boundaries of the three regions can be easily calculated using the lens equation ([6]) by 
imposing 9± = Othr = 3 mas in the case of GRAVITY and solving for 7 as a function of -Dls- I n 
Figure [2] we can see the shapes of the three regions in the plane (7, -Dls)- 

Region I covers the immediate neighborhood of Sgr A*, when the source is too close to the 
black hole to generate a distinguishable primary image. Using the weak field/thin lens formalism 
(Eq. flOl) . an approximate expression for its boundary is 

-DOL + ^LS/) 2R S . , 

7/ -" = ~ ^ ~ D^Mr (16) 

The tip of region I at 7 = extends up to 



3.3. Three regimes where to look for lensing in the Galactic center 




(17) 
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Beyond this point the primary image will always be resolved by GRAVITY whatever the 
alignment angle 7. The weak-deflection analytical expression tracks the exact boundary extremely 
accurately, as # t hr — 3006*5. 

The possibility to observe both images as separate objects from Sgr A* is restricted to sources 
within a cone starting at the tip of region I. The analytical expression coming from the weak 
deflection approximation is just the opposite of Eq. (|16p 



D l + D ls 2R s 

111-111 = n U thr + — . (18) 

-^LS ^OLUthr 

For sources at large distances we have 7/1-77/ — »■ 0.39°. 

In the next three sections we will closely examine the three regions just defined, pointing to 
possible interesting observational perspectives with GRAVITY, considering either photometric or 
astrometric detection channels. 



4. Region I 

In region I of the parameter space (Figure [2]) , the primary image of a generic stellar source 
is within 3 mas from the black hole. Then, it cannot be resolved from the radiation emitted by 
Sgr A* and forms a single blend along with the secondary image. Obviously, in this situation it is 
difficult to extract pure gravitational lensing signals on the lensed star, as everything is mixed with 
the variable flux from Sgr A*. 



4.1. Astrometry 

Let us consider astrometric variations first. The centroid of the blend composed of Sgr A*, 
primary image and secondary image of the lensed star can be calculated by evaluating the position 
of the images by Eq. ([6]) and their magnification from Eq. ([7]). Of course, the weight of the two 
images with respect to Sgr A* depends on the relative fluxes in the if -band. By considering that a 



magn itude of K = corresponds to 667 Jy and that the extinction is = 2.8 (jEisenhauer et al 



2005), Sgr A* has an apparent magnitude of K ~ 17 in its quiescent state, corresponding to 1.1 
mJy As mentioned before, this flux is highly variable even within the same day. However, for 
the sake of simplicity, we have taken -KsgrA* = 17 for the calculation of the centroid shift. As for 
the stellar source, typical values of S-stars range from K = 13 to K = 17. This also means that 
if gravitational lensing effects take place, the centroid will be typically dominated by the primary 
image shift, as Sgr A* is relatively faint. 

In Figure [3] we show the expected astrometric shift of the centroid in region I due to gravita- 
tional lensing effects for a source with K = 15. In practice, we are comparing the centroid position 
of the blend including Sgr A* and the two lensed images, with the centroid position obtained 
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with Sgr A* and the unlensed source. As expected, the shift is higher for 7 — > 0, i.e. for better 
alignments. Considering that the astrometric sensitivity of GRAVITY is of the order of 30 //as, 
gravitational lensing effects would be detectable for a large portion of this region up to 7 = 30°. 
This anticipates the results of the astrometric effects in region II, to be discussed in the next section, 
where the primary image is resolved and things are much cleaner. In region I, instead, the variabil- 
ity of Sgr A* would add some noise on the astrometric measurements. A reliable extraction of the 
gravitational lensing astrometric shift would require follow-up observations of the whole transit of 
the source behind Sgr A*, so that the noise due to Sgr A* can be efficiently subtracted. 

Strong field and spin effects on the total centroid shift in region I are tiny, staying below the 
/ias even for sources closer than 10 pc. 

4.2. Photometry 

The second channel for the detection of gravitational lensing effects comes from photometry. 
In this case as well, the variability of Sgr A* adds noise to the flux variations of the primary and 
secondary images. We have thus fixed a threshold for photometric detection of lensing magnification 
effects of AK = 0.1. The time-scale of lensing variations depends on the radius and the eccentricity 
of the orbit of the source S-star. For example, the transit of a star on a circular orbit at a distance 
Z>ls = 0.01 pc behind the black hole in the region within 3 mas of Sgr A* lasts about 100 days, 
scaling with v^Dls- As this time is much larger than the scale of the variability of Sgr A*, the two 
phenomena can be disentangled with suitable follow-up observations. 

Figure H] shows region I split in two zones: the left one where the photometric variation due 
to gravitational lensing is detectable and the right one where the magnification effect is too small. 
The three curves track this boundary for three values of the source magnitude in JC-band (13, 15, 
or 17). Note that an S-star with 7 > 10° will not induce any perceptible flux variation even if it 
were more brighter than S2. Based on our knowledge of the orbits of several S-stars, we can say 
that GRAVITY will always be able to separate the primary image from Sgr A*. This means that 
at present we do not know S-stars falling in region I. However, GRAVITY will likely discover more 
stars that are too close to be detected with present facilities. It appears a reasonable guess that 
some of them will enter the regime described in this section along their orbit. If this is the case, 
their gravitational lensing effects can be studied along the lines traced in this section. 

Interestingly, the photometric magnification of a stellar source provides a first channel for the 
detection of the spin of the black hole. In Fig. [J] the dotted line bounds the region in which the 
displacement of the Kerr caustic discussed in Section 13.21 modifies the magnification in such a way 
that the total observable flux changes by more than AK = 0.1. Of course, in order to distinguish 
the spin signal, we need to know the relative positions of Sgr A* and the source S-star to very high 
accuracy from previous observations. Contrarily to the spin, Post-Newtonian terms in spherically 
symmetric black holes do not give rise to detectable effects. 
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We finally note that the threshold precision of 10% that we have requested here has been chosen 
by considering the scatter in the data collected with the Naos-Co nica (NACO) system m ounted at 
the fourth unit telescope Yepun of the VLT (see over Figure 8 of iGillessen et al.1 l2009al fl We are 
quite confident that GRAVITY should be able to achieve a similar precision in measurements of 
the same type. 



5. Region II 

In this region, the primary image of the source star is detached from the black hole, whereas 
the secondary image is not. Obviously, this regime includes the no-lensing limit with j3 ^> 6e, in 
which the primary image is basically unlensed and the secondary image becomes negligible. All 
known S-stars fall in this region, so that GRAVITY will be able to resolve their primary image all 
along their orbits (in current VLT observations, instead, S2 and S6 used to be unresolved from Sgr 
A* for some periods). 

In this regime, we can look for four different gravitational lensing effects: astrometric shift of 
the primary image, magnification of the primary image, centroid shift of the blend composed of 
Sgr A* and the secondary image, magnification effects on the same blend. 



5.1. Astrometry on the primary image 

The astrometric shift on the primary image due to gravitat i onal l ensing by the central black 



hole has already been studied in detail by iNusser &: Broadhurst] (120041 ) . These authors could not 



envisage the breakthrough in astrometric accuracy promised by GRAVITY, so that their conclusions 
were not very optimistic for the observational perspectives of such an effect in a short term. 

For this reason, it is interesting to revisit the expected astrometric shift in the source coordinate 
space in the light of GRAVITY expectations. Figure [5] shows contour lines for the astrometric shift 
at several values. Strikingly, many of the known S-stars would have experienced a measurable 
astrometric shift during their passage behind Sgr A*. Even at alignment angles as large as 7 = 35° 
the astrometric shift is of the order of 30 //as. This encourages to look for astrometric shift effects 
on known S-stars as soon as GRAVITY becomes fully operational. To this purpose, in Figure [6] 
we show the position of the S-stars in the (7, .Dls) plane in the years from 2014 to 2024. We see 



2 The question of how accurate the K values are, is actually difficult, because is not easy to understand how 
the systematic uncertainties affect the measurements. For instance, the extinction around S gr A* is not constant 
spatially but at the contrary, as showed by the accurate extinction map of ISchodel etafl (|2007l l. is very patchy. This 
means that stars that hang around the Galactic center experience higher and lower extinctions, which modulates the 
measured fluxes. One of the possible explanation of the S2 flux increasing observed in 2002, when the star passed its 
pericenter, is just that (Gillessen 2009, private communication). 
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that S17 will have an appreciable shift for a few years, while several S-stars may have a marginally 
detectable shift. Interesting targets will be S29 and S2, which will make rapid incursions in this 
plot as they will pass through their periapse. 

Figs. [5] and [6] clearly show that the shift is practically independent of the source distance. 
This statement can actually be supported by a simple analytical derivation. In fact, since light 
rays of the primary image pass very far from the event horizon, the weak field lens equation ([8]) 
gives the same results as the exact lens equation ([6]) to very high accuracy (the difference at most 
reaches 5 //as for 7 — > 0). Conversely, at large alignment angles 7, the small angle approximation 
(jlOP fails. Starting from Eq. (|8j), solving perturbatively for u in powers of Rs and taking the limit 
D h > D LS , we find 

A>l sin 7/2' v ; 

which is independent of Dls- 

Proceeding in a similar way from the small angles thin lens equation (|10p . we get 

A6 ~ ^L, (20) 

which coincides with the small 7 expansion of Eq. (|19p . but gives an error of order 10% at 7 = 30° 
already and thus cannot be used for accurate shift calculations on S-stars. 

We finally note that this gravitational lensing shift should be taken into account in any ac- 
curate orbit determination undertaken from GRAVITY data. This effect will help setting further 
constraints on the black hole mass, the mass distribution around it and the orbital parameters of 
S-stars. As anticipated, Post-Newtonian and spin effects for the primary image shift in region II 
are negligible. 



5.2. Astrometry on the secondary image 

In region II the secondary image is still blended with the emission from Sgr A*. As the source 
aligns with the black hole, the secondary image gets brighter and moves away from Sgr A*. Both 
these effects contribute to shift the centroid of the blend in the direction opposite to the primary 
image. Obviously, the centroid shift is determined by the balance between the luminosity of Sgr 
A* and that of the source S-star. 

Similarly to what we did for region I, we calculate the centroid shift of this blend. In Figure 
[7]we show the detection limits of GRAVITY (set to 25 //as) for three different values of the source 
magnitude in the iT-band. Interestingly, S27 and S6 at their closest alignment would have fallen 
inside the detection zone if they were more luminous (S6 has K = 15.4 and S27 has K = 15.6). No 
other known S-stars generate secondary images prominent enough to be detected by GRAVITY. 
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Post-Newtonian effects on the secondary image are larger since this image is formed by light 
rays passing close to the event horizon. However, in the region of interest in Figure the difference 
between the exact centroid shift and the shift calculated by the weak field formalism stays below 1 
//as. The spin effects reach at most 2 //as, one order of magnitude below GRAVITY'S sensitivity. 
In principle, using a very large sample of observations of different stars whose orbit is very well 
known, it might be possible to estimate the spin and Post-Newtonian terms by a careful statistical 
analysis. So, this is an interesting possibility opened by a measure that is in principle feasible 
with GRAVITY, provided that some S-stars pass through this region. Obviously, the lensing signal 
is confused here with the variability of Sgr A*, which needs to be subtracted with long-period 
observations. 

5.3. Photometry on the primary image 

The magnification of the primary image is another possible channel for detecting gravitational 
lensing. It is inversely proportional to the alignment angle (as long as the source radius is negligible), 
so that it requires quite accurate alignments. The expected magnification for the known S-stars is 
very low (we estimated that the best case is // = 1.01 for S6) and so surely undetectable also by 
GRAVITY. If we consider a detection threshold of AK = 0.1, the candidate sources must be on the 
left of the dashed line in Figure El i.e. within 7 = 1°. It is a small window for this channel, but not 
impossible to achieve for S-stars to be discovered yet. By crossing astrometric and magnification 
measurements it is possible to get even more stringent constraints on the lens nature. 

5.4. Photometry on the secondary image 

The photometry on the secondary image is complicated by the blending with Sgr A*. The 
threshold at AK = 0.1 allows for a relatively wide detection region depending on the source 
magnitude. In Figure [8] we see its boundary for three values of the magnitude of the S-star. With 
respect to astrometry, the region extends at larger angles at small -DlSi while it becomes thinner at 
larger Z?ls- The difference between the exact and weak field values is smaller than in astrometry 
and is thus not shown in the figure. Also in this case, no known S-stars fall within the detection 
limits of GRAVITY. 

Note that the detection limit for K = 17 practically coincides with the limit for the primary 
image. This is a trivial consequence of the fact that the magnitude of Sgr A* in its quiescent state 
has been taken to be K = 17, i.e. equal to the unlensed source magnitude. 

Post-Newtonian and spin contributions to photometric amplification remain small both for the 
primary and for the secondary image, starting to be significant only for 7 < 1'. 
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6. Region III 

In this region, the two images of the source star appear distinct from Sgr A*. Obviously, this 
is the most favorable case for a gravitational lensing detailed study. The mass distribution and 
the distances can be constrained by simultaneous measurements on the primary and the secondary 
image position and magnification, so that the greatest benefits for a study of the central environment 
is achieved. Unfortunately, region III is a cone with aperture angle of 7 = 0.39° only, so that the 
chances for such a good alignment are quite low. 

Furthermore, we still have to discuss the actual detectability of the secondary image, which 
is always much fainter than the primary. Following our discussion of GRAVITY performance 
in Section 12. 4| we consider a threshold magnitude of K = 19 for the secondary to be clearly 
detected. Depending on the source intrinsic luminosity, in Figure [9] we draw three detection limits 
corresponding to this threshold. These limits demand an even higher degree of alignment for the 
secondary image to be effectively visible. 

In the weak field/thin lens formalism, these curves can also be calculated analytically by im- 
posing that the magnification of the secondary image (|13p is above the threshold /i = io~°- 4 ( 19- ^) , 
where K is the magnitude of the source. We obtain 



/ 2fl s (l + 2/x)-2V/x(l + /i) 

Consider that at a given distance -Dls> the radius of the detection zone scales as -DlsT — V^LS) 
i.e. the detection zone effectively grows with the distance, contrarily to the impression given by 
Figure [H 

Another interesting aspect is that throughout region III the contribution of Post-Newtonian 
and spin effects to the astrometric shift is relatively high: 5 //as for the primary image and up 
to 10 //as for the secondary image. These values are still below GRAVITY sensitivity. However, 
as already noted before, a large number of observations can be used to study these effects with a 
statistical approach. 

We finally remind the estimates of previous works on the number of bulge stars that should be 
lensed by Sgr A* at a given time with a threshold of K = 23 for the sec ondary image is around 10 



()Jaroszynskilll998al : lAlexander fc Sternbergill999i : IChaname et al.ll200ll ). Here we consider a more 



conservative threshold for GRAVITY, since we are considering its specific use for interferometry. 



7. Discussion and conclusions 



Thanks to the coherent combination of the four 8-m Unit Telescopes and the four 1.8-m Auxil- 
iary Telescopes, which can be relocated on 30 different stations, the VLTI is the only interferometer 
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that allows to perform very high-precision photometry, with mas angular resolution, and a total 
collecting area of 200 m 2 . In particular, the size of the total field of view combining the four 8-m 
Units is of 2". Since no other interferometer of such class are currently in the desi gn phase, the 



VLT I will be at t he leading-edge sti ll for many years to come. After the VI NCI. MIDI (ILeinert et al 



2011 



20031 ). AMBER (jPetrov et al.ll2007l ) and the visitor instrument PIONIER (|Le Bouquin et al. 
GRAVITY arises as the second generation VLTI instrument and will perform high precision narrow 
angle astrometry as well as interferometric imaging in the i^-band (2.2 /mi). 

Taking into account the capabilities of GRAVITY, we have analyzed the possibility to observe 
gravitational lensing effected generated by the massive black hole, located in the geometrical center 
of the Milky Way, on the S-stars orbiting around it. 

We have divided the parameter space into three regions (Figure [2]), depending on the resolution 
of the primary and the secondary image from the black hole and we have carefully investigated 
astrometric and photometric channels for the detection of gravitational lensing effects. 



Region I For stellar sources in the immediate neighborhood of Sgr A*, both the primary and the 
secondary images are confused with Sgr A*. 

Gravitational lensing is responsible of a shift in the centroid with respect to the unlensed 
configuration. This shift is measurable by GRAVITY up to an alignment angle of 7 ~ 30° 
(cfr. Figure [3l see Figure [1] for the geometrical definition of 7). 

Photometric detection of a magnification effect is possible in the inner region 7 < 2-^5° 
(Figure H} . 

These effects have timescales much longer than the daily variability of Sgr A*, which could 
be subtracted by a dense sampling. 

Region II If the S-star is far enough from the black hole, the primary image is resolved by GRAVITY. 
All known S-stars fall in this region. 

In this case, the astrometric shift of the primary is appreciable by GRAVITY for many known 
S-stars (Figs. [MB]) in a cone with aperture 7 ~ 40°. This shift must be taken into account in 
the computation of orbital parameters and may help constraining the mass of the black hole 
and of its environment. The star SI 7, in particular, will have the highest astrometric shift 
for several years. 

The magnification effect on the primary is more difficult to see and requires 7 < 1° (Figure 
E]). 

The secondary image is blended with Sgr A*, but is able to yield a measurable centroid shift 
in the zone between 1° < 7 < 8° depending on the source distance and magnitude (Figure [7]). 

The magnification effect on the secondary image might be easier to detect (Figure [8]), given 
the relatively low luminosity of Sgr A* in its quiescent state, but it is affected by the problem 
of the variability of the comntaminant emission, as mentioned before. 
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Region III In a cone of aperture 7 = 0.39°, the stars will generate two images well separated from Sgr A*. 

Of course this is the best situation for investigations using gravitational lensing to determine 
masses and distances in the Sgr A* environment. By imposing that the secondary image is 
not too faint, we further restrict the zone suitable for such studies (Figure [9|). 

For sources exceptionally aligned on the optical axis and close enough to the black hole, the 
magnification is sensitive to the spin of the black hole (Figure H|). Apart from this case, spin 
and Post-Newtonian corrections to the weak field lensing are generally below the sensitivity of 
GRAVITY in all examined regions. However, the corrections to the astrometric shifts in Region 
III and the shift of the secondary image in Region II are not far from GRAVITY sensitivity and 
may be extracted through a careful statistical analysis on a large number of observations. 

In conclusion, we can say that GRAVITY will be able to provide the first evidences of gravi- 
tational lensing effects by Sgr A*, thanks to its striking astrometric abilities. The astrometric shift 
due to gravitational lensing must be taken into account in the reconstruction of the orbits of S-stars 
and can provide independent precious information on the physics of the massive black hole at the 
center of our Galaxy. New exciting windows on effects beyond the weak field gravity might also be 
opened if more S-stars are discovered. 

We wish to thank Stefan Gillessen for useful discussion about the photometric accuracy of NIR 
measurements. 
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